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Abstract 
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I. INTRODUCTION 



The production of an electron-positron (e + e~) pair by a photon in a strong atomic field 
has been investigated already since many years both theoretically and experimentally be- 
cause of the importance of this process for various applications, see Refs. [l|, |2J. The cross 
section of this process in the Born approximation is known for arbitrary energy u of the in- 
coming photon, Refs. {3, ^] (we set h — c — 1 throughout the paper). The effect of screening 



in this approximation can be easily taken into account using the atomic form factor [5|]. For 
heavy atoms, it is however necessary to take into account the Coulomb corrections. These 
corrections are higher-order terms of the perturbation theory with respect to the parameter 
Za, where Z is the atomic charge number and a = e 2 ~ 1/137 is the fine-structure con- 
stant, with e being the absolute value of the electron charge. The formal expression of the 
Coulomb corrections, exact in Za and uj, was derived in Ref. This expression has a 
very complicated form which leads to difficulties in numerical computations. The difficulties 
grow as oj increases, so that numerical results have been so far obtained only for u> < 12.5 
MeV, Ref. Q- 

In the high-energy region u ^> m (m is the electron mass) , considerations become greatly 
simplified. As a result, a simple form of the Coulomb corrections was obtained in Refs. js), Q 
in the leading approximation with respect to m/u. However, this result has good accuracy 
only at energies u> > 100 MeV. The theoretical description of the Coulomb corrections for 
the total cross section at intermediate photon energies (5 -j- 100 MeV) was based during 
a long time on the "bridging" expression derived in 10] . This expression is actually an 
extrapolation of the results obtained at u < 5 MeV. Results for the spectrum of one of the 
created particles at intermediate u were practically absent. Recently, an important step was 
made in Ref. 11( where the first corrections of the order of m/u to the spectrum as well 



as to the total cross section of e + e~ photoproduction in a strong atomic field were derived. 
The correction to spectrum was obtained in the region where both produced particles are 
relativistic. It turns out that this correction is antisymmetric with respect to replacement 
e + -H- e_, where e + and e_ are the energy of the positron and the electron, respectively. 
Since the correction to the total cross section resulted to be very large, it is not related to 
the central region, where the created electron and positron are ultrarelativistic, and it comes 
from the region close to the end of spectrum where e + ~ m or e_ ~ m. In Ref. 11], the 
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correction to the total cross section was obtained with the use of the dispersion relation for 
the forward Delbriick scattering amplitude but not by the direct integration of the spectrum. 
Note that the account for the correction to the total cross section leads to good description 



12|. 



of available experimental data at intermediate photon energies, Refs. 

In the present paper, we calculate the electron (positron) spectrum in the process of 
e + e~ photoproduction in a strong Coulomb field in the case e_ ~ m (e+ ~ m) and u 3> m. 
We show that the Coulomb corrections drastically differ form that obtained in the region 
where e + 3> m and e_ 3> m. In an analogous way, we have also derived the spectrum of 
bremsstrahlung in a strong Coulomb field in the region where the radiated photon has the 
energy close to that of the initial electron. 



II. GENERAL DISCUSSION 



Following the usual Feynman rules (see, e. g., Ref. [13]), the e + e~ photoproduction 
cross section, at leading order in the interaction between the photon field and the electron- 
positron field, averaged over the polarization of the incoming photon and summed up over 
polarizations of electron and positron has the form 



M" = J dxU p _, x _(x)rV p+M (x)e^ x , (1) 

where (00, k) is the four-momentum of the photon, (e_,p_) and (e+,p+) are the four- 
momenta of the electron and the positron, respectively, and A_ and A + are their polarization 
indexes. Also, U p _ ) \_ and V p+: \ + are the corresponding positive-energy and negative-energy 
wave functions in a strong Coulomb field and 7 M are the Dirac matrices. In the following, we 
will calculate the spectrum, i.e., the cross section integrated over the angles of the vectors 
p + and p . Due to rotational symmetry, this quantity is independent of the direction of the 
photon momentum k. Therefore, we can average it over this direction, i.e., integrate both 
sides of Eq. (P) over dOfe/(47r). This results in the replacement 
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where R = \y — x\. We can also use the known relations \l£ 

j^U r , x (x)U r , x {y) = ^-5G(x,y\ e ), 

where e > and 5G(a;, j/ 1 ±e) is the discontinuity on the cut of the electron Green's function 
in the Coulomb field, to obtain 

" 1 f dxdySp[^SG{x 1 y\€.-u)^G{y 1 x\€-)]^^- 1 (4) 



de- 2u J J ' " uR 

Since the representation of the quantity SG(y, x |e_) is independent of the basis em- 
ployed, we can also write 5G(y, x |e_) in the form 

8G(y, x |e_) = U j>Vtlt (p-, y)U j)fftlt (p-, x) , (5) 

and represent the spectrum as follows 

le~ = 2^~^Z / / dx d yUj,«AP-> x)Y5G(x, y |e_ - u)^ p U j>(Tjlx (p^, v) ^^^~ » ( 6 ) 

where /3_ = p_/e_, a = ±1. Here, Uj^^(p, x) is the positive-energy wave function with total 
angular momentum j, parity (— l)-^ "/ 2 , and projection /x of the total angular momentum 
along some quantization axis. The explicit form of this function is presented in Appendix 

m 

Below, we assume that u 3> m. If the momentum of the electron is p- ~ m, then 
the formation length of the process is of the order of the Compton wavelength 1/m and 
the positron is ultrarelativistic with a typical angular momentum l + ~ oj/m ^> 1. This 
circumstance allows us to use the quasiclassical Green's function obtained in Refs. [141 

n n 

15| starting from a convenient integral representation derived in Ref. [16] of the exact 



Green's function of the Dirac equation in a Coulomb field. For the reader's convenience, 
we present the formula for the discontinuity of this Green's function in Appendix |AJ Eqs. 
( IA8I) . Moreover, it can be seen that if p_ <C p+, then the main contribution to the integral 
in Eq. ([6]) is given by the region where the angles between vectors x and y are not small 
(or close to n). This is the region which also gives the main contribution to the bound-free 



photoproduction cross section at oj ^> m, see Ref. 



17j . In this region the expression in Eq. 
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(1A8I) becomes essentially simpler, see Eq. f ]A9j) . By substituting Eq. (1A9[) in Eq. ([6]) and 
by keeping the terms which do not contain highly oscillating functions, we arrive at 
da a v-^ f f dxdy- 

dT = AncoB ^ J J ~tf~ U i^' p -' X >( 1 cos< ^ ~ 11 ' nsm W>,/AP-, V) > ( 7 ) 

where = e_i? + 2Zas (the notation is explained in Eq. (1A9|) ). By employing the explicit 
form (1A1I) of the electron wave function, we finally arrive at the following expression of the 
cross section of pair production in the case of slow electron 

da a 



de- 4ir 2 u8. 



^0 + 5)//Sf (Fcos *- rsm ' 



F = f(x)f{y)P l {t) + g(x)g{y)Mt), 

T = ^{f{x)g{y)[xP v {t) -yPi(t)\ + g(x)f(y)[yP ll (t) - xP t {t)]} , (8) 

where t = x ■ y/(xy), I = j + a/2, I' = j — a/2, and Pi(t) are the Legendre polynomials. 
The cross section on the other end of the spectrum, i.e. in the case of slow positron with 
momentum p + ~ m, is given by Eq. (jSJ) with the replacement e_ — » e + , f3_ — > j3+ = p+/e + , 
and Z — > —Z. With the same procedure, one can also derive the spectrum of bremsstrahlung 
by a high-energy electron for the case where the final electron with momentum pi and energy 
e-y is slow (pi ~ m). It turns out that this spectrum is given by the same formula (|HJ) with the 
obvious substitutions p_ — » p\ and e_ — > e±. Note that the correction to the bremsstrahlung 
spectrum of the order of m/e\ in the case of initial and final electrons both ultrarelativistic 
was obtained recently in Ref. HI] from the corresponding results for pair production, and 
in Ref. 18| directly from the matrix element of bremsstrahlung. 

It can be shown that three integrations in Eq. ([8]) can be performed analytically and the 
spectrum da/de_ becomes 

i 



da 2a 



de „ zZO + l) j j dxd y [ dt^{Fcos(j>-Tsm<j>). (9) 

J> -1 

It is convenient to multiply the integrand in this formula by unity written in the form 

i 

1 = / du2u5 (u 2 - - - ) , (10) 







{x + y) 



to change the order of integration over the variables t and u, and to take the integral over t 
(by exploiting the 5-function). After that, we pass from the variables x and y to the variables 



p and v such that x = p(l + v)/2 and y = p(l — v)/2. As a result we obtain 

oo 11 

tt = ^kv{ i+ k) I pdp S dv S v' Fcos *- rsin *) 

• J> -1 \v\ 

F = f(x)f(y)P l (t )+g(x)g(y)P ir (t ), 
T = ^{f(x)g(y)[(l + v)P l/ (t )-(l-v)P l (t )] 

+g(x)f{y)[(l - v)P v (to) - (1 + v^ito)} 
1 + v 1 — v 1 — u 2 

X = p^-, y = p—^-, t = 2^—^-1, 

$ = e_wp + Za In ( ■ (H) 

This formula is still not convenient for numerical calculations because of the strong oscilla- 
tions of the integrand in the vicinity of the point u — 0. In order to overcome this difficulty, 
we write 

F cos $ — T sin $ = ReM(w, v, p) , M(u, u, p) = (F + iT)j* . (12) 

Then, by using the properties of the integrand M(u,v, p), we make the following transfor- 
mation 

1 1 oo 1 1 oo 

dv J — J pdp M(u,v, p) = — J du J dv J pdpM(u,vu,p) 

1 \ v \ -1-10 

oo 1 oo 

i 
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Jdrjdvjpdp M^ — l + ir,v(—l + ir),p) - M (l + ir,v(l + ir), p^j . (13) 
o-io 

In the first step we changed the integration order of the variables u and v and then we 
performed the change of variable v — > vu. In the second step, we changed the contour of 
integration, by exploiting the fact that the contribution along the straight path from the 
point u\ = 1 + ioo to the point 112 = — 1 + ioc vanishes due to the exponential function 
exp(i$) (see Eq. flTTJ and also Fig. 1). This form of integral is appropriate for numerical 
calculations. Note that the integral (TT3"j) has zero imaginary part, so that it is not necessary 
to take real part of it afterwards. 
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III. ASYMPTOTICS /3_ -> 



Let us consider the spectrum Eq. (TTTT) in the limit /?_ = p-/e_ — >■ 0. Substituting the 
asymptotics ( 1A5j) of the functions /(r) and <?(r) in Eq. (jlip . we obtain 

1 i 



da air 



de- oom 2 (Za 



+ V / ^ y ^(^ocos$ -T sin$ ) , 



F = {[k 2 + (Za) 2 ] J^( m ) J 2 ,( V2 ) + ^4 y (vi)4 7 (V2)}[P j+ i/2(to) + Pi-i/2{to)\ 

~\ (j + ^ [^2 7 (^l)^2^2 7 (^2) + ^27(^2)^1^7(^1)] [Pj+l/2(to) ~ Pj-l/2(t )} , 

To = ^{v[J 2l {rii)ri 2 J' 2l {ri 2 ) - J 27 M^i4 7 (^)][ P i+i/2(*o) + Pj-i/ 2 (to)] 

-4 (j + ^ J2 7 (^l)J2 7 (^)[Pi+l/ 2 (to) " ^-l/2(to)]} , 

771 = 2a/Z«p(1 +w) , r/ 2 = 2a/ Zap(l - u) , 

/l + u\ 1 — u 2 

% = up + Zalni- , t = 2- r-1. (14) 

VI — u/ 1 — v A 



The integration over the variable p can be taken by using the relation 19] 

00 22 
J dxe^J^a^J^by/x) = exp - i - ^ J , (15) 



The remaining integrations over the variables v and u can be performed numerically by 
employing the transformation (fl3|) . The largest contribution to the sum over j is given by 
the term with j = 1/2. The contribution of the term with j = 3/2 is essentially smaller, while 
that with j = 5/2 is less than one percent even for large Z. In our numerical calculations 
here we have included terms with j = 1/2, j = 3/2 and j = 5/2. The results for uda/de- 
in units of a = a(Za) 3 /m 2 at zero electron velocity is shown in Fig. [2] as a solid curve. At 
Za — > we obtain uda/de- = Air a in agreement with previous results. In Ref. 20) by Deck 
et al., the following formula for uda/de- at zero electron velocity was suggested 

da^ aiZaf 2ixZa ( 4tt „ \ 
u— = 4tt^— ± — — — — 1 - — Za . (16) 



de_ m? exp(27rZo;) — 1 \ 15 

This formula is shown in Fig. [2] as a dashed curve and it is clear from the figure that Eq. ( 1T6"|) 
is applicable only at small values of Za. Also, note that the cross section udas/de^ in the 
Born approximation vanishes in the limit /3_ — > 0, since at /3_ 1 it scales as udas/dt- ~ 
2a(Za) 2 (3-/m 2 . In Ref. [7| the results for cuda/de- was obtained at u — 40MeV and 



e_ = 1.008 m (which corresponds to /3_ = 0.1265). These results are shown in Fig. [2] as 
a dotted curve starting from Z — 11. One sees an excellent agreement of our results with 

n 

those of Ref. (7J. At Z = 1 there is disagreement because at small Z and /?_ = 0.1265 the 
contribution of the Born term is also important. Finally, we observe that, as expected, the 
spectrum of positron at small positron velocity tends to zero because in this case the wave 
functions are exponentially small (see Appendix [A]) . 



IV. CROSS SECTION AT NON-ZERO ELECTRON VELOCITY 



In order to obtain the spectrum for non-zero electron velocity, we substitute the explicit 
form of wave function ( lAlj) in Eq. (fTTj) and use the relation F(a, (3-,x) = e x F(/3- — 
a, /3_, —x) for the confluent hypergeometric function. We come to the following expression 
for the cross section at p_ ~ m 



da 
M -- 



a 



L=l 



|r(7+ 1 + iv 

[r(2 7 + 



^^Re J du J dv J dpp^ +1 (l-v 2 u 2 ye^M 



F\F 2 
7 — iv 



m 

iv—z-A. 

61 



L 1 + 



a 



+ 



FiF 2 (l + /9_t;) + FiF 2 (l - P-v) 



A 



+ ) 



F\F 2 



iv^ r A + + L 1- — 

61 \ U 



Fi,2 = F{j - iu, 2 7 + 1, -ip(l ± vu)) , F 1:2 — F(j + 1 — iv, 2 7 + 1, -ip(l ± to)) , 
A ± = P L (t)±P L _ 1 (t), $=(^ + l)p + Zaln(l±^) , t = 2^-^-1, (17) 

where v = Za/ (5_ and 7 = \J L 2 — (Za) 2 . Then, we perform the transformation ( TT3l) 
and take analytically the integral over the variable p using the relation (see mathematical 



Appendix f in Ref. 



2l|) 



OC 



r( 7 )A a+a '- 7 (A - k)~ a (X - k')- a 'F ( a, a', 7 



kk! 



where F(a, b, c, x) is the hypergeometric function. After that we take numerically the inte- 
grals over the variables v and u. 

The cross section in Eq. ( ITT)) can be represented as 



da (icr £ 



L=l 



(19) 



Unfortunately, the convergence of the series (fT9|) is not fast, and it is necessary to take into 
account terms do^jde^ with rather large L, especially at e_ ^> m. In order to overcome this 
problem, we make the following approximation 

^dc L ^ (d<r L da£\ A A _^da A 

where Lq is a large integer and da A /de- is given by the expression for dai/de- with the 
replacement 7 = >J L 2 — (Za) 2 — » L. The convergence of the series ^lLiI^ "^/^ 6 - ~~ 
da A /deJ) is much faster than the convergence of X^lLi da^/de^^ and it is not necessary to 
take very large value of L Q (we have seen that by choosing L = 5, an accuracy of about 
10 % is reached), while the series T> A can be summed analytically Since the summation is 
not straightforward, we report some steps in the next paragraph. 



A. Calculation of T, A 



It is convenient to perform calculation of S A starting from Eq. fl4]). After the replacement 
7 — > L, the expression for the discontinuity of the electron Green's function is given by Eq. 
( 1A8I) . By using this expression and by passing to the same variables as in the derivation of 
Eq. (jlip . we can again employ the convenient transformation fll3p . In this way, the integral 
over the variable v can be easily performed and one obtains 



1 a 



du 



M A {-1 + iu) - M A {1 + iu) 



M A (u] 



3 A>$> 



dp p e 



dr 



sinh 2 t 



x 



exp [i(2z/r + p cothr)] 



1 / u \ 

— 1 J (w) - 1%Za(\ - u 2 ) cothr 

p_ V 3 / w 



2u , T . . . pu(l — u 2 ) Ji(w) 
+— cothrJ H- r \ ; u J 
o 3 sinh r w 



w 



p\J\ — M 2 



$ = ^- + Za In 



(21) 



sinh r /3_ \_ 1 — u / 

Since the contour of integration over r passes in the positive direction around the point 
r = 0, we can make a shift r — > r — in/ 2 and take the integral over r in M A {u) using the 
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relation 



22] 



dr 



coshr 



exp (2Ar — b tanh r) J; 



2n 



c - fe3 2M r (l + A + / i ) r (|- A + / i ) 



coshr; [T(2//+l)] 2 

— oo 

xF U - A + i 2// + 1, 6+ V& 2 - a 2 ^ fL-A+^,2/1 + 1,6- V& 2 - a 2 

Then we take the integral over p with the help of Eq. (fl8l) . Finally, we have 

,2\ 2 



M (m) = 2i 



1 



+ Za(l — u ) uv 

3 



r (2) 
'11 



2i 



u 



u 



+(l-w)t-(l-iv) I 1-- \lg>--(2-w) I 1- - I [Jii'+J- 



22 



'It 



r(l) , r(l) 



+ | (2 _ f , )(4 ;»_ / o.) + i( 1 -^)( 



iff + iff) - 



31 ~ J 13 

2 



1 



ir 



+ 



2m 





1 


+ 






~2Za 


+ 






u 


® = ( 





\ 3 , 

(l-u 2 )+ % -Liv+ ^-5 



r(2) _ r(2) 
21 i 12 



l-(l) _ l-(l) 
J 21 J 12 



"I) 



" A A=-t(l-Ht/0-) 

I jk (\) = \i+ k - 2i »(\ + i(i + u )) iu - j (X + z(l - m))^ 
xF I j — k — iv, 2 



(22) 



r (2) , r(2) 
J 12 "r -<21 



(A + i(l-«))(A + i(l + w)) 
This expression is particularly suitable for numerical integration. In the next section we 
report our results obtained starting from Eqs. (1171) in the approximation ( 1201) . 



(23) 



V. RESULTS AND DISCUSSION 



The cross section dos / de_ of the pair production process in the Born approximation is 
well known (see, for example, Ref. Il3|). In our limit (u ^> m and p_ Cu), it has the form 



das c"o 2e_ 



2e_p_ In 



m 



2 2 1 2 / e - + P- 
p_ - m In 

m 



^0 = — -tt— • (24) 



In particular, at p_ <C m it is das/de^ = 2a p-/m, while at p_ ^> m the same cross section 



has the asymptotics 



do B _ 4^0 



In 



2e_ 



m 



(25) 
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Now, the leading Coulomb correction to da/de- at m <C e_ <C u reads 



191 



da, 



(0) 

c 



4a 



/(Za), f(Za) = Re[i/;(l + iZa) + C}, (26) 
= 0.577... is the Euler constant. The correction (1261) is 



where ip(x) = dlnT(x)/dx, C 
independent of e_ and it is the same for electron and positron (i. e., it is an even function of 
Za). The next-to- leading correction was calculated recently in Ref. and at m <C e_ <C uj 
it has the form 



da. 



(i) 
c 



a 7i m 



Reg(Za), g(Za) = Za 



T(l-iZa)T(l/2 + iZa) 



(27) 



de_ w 2e_'"^ v " ^ v ' r(l + iZa)T(l/2 - iZa) 

This correction has opposite sign for electron and positron since g(x) is an odd function of 

Za and it increases the cross section for slow electron while it decreases it for slow positron. 

In Fig. 3 and Fig. 4 we show our results for the Coulomb corrections toa^dac / 'de T = 

uaQ 1 (da / de T — daB/de T ) to the spectrum for slow electron and positron, respectively at 

different values of Z (continuous curves). These results are compared with the asymptotic 

expressions ua 1 da^ /de T (dashed curves) and uaQ l (da { § > /de T + da { ^ > /de T ) (dotted curves). 

On the one hand, one can see that for each Z our results tend at large energies to the constant 

value —Af(Za). On the other hand, the next-to-leading correction da^ /de T essentially 

improves the agreement between exact results and asymptotic ones both for slow electron 

and positron. 

In Ref. 11] the next-to leading correction to the total cross section, a^ oc m/u, was 



also obtained. It reads 



a. 



(i) _ ma 



c 



UJ 



7T 



-—luig(Za) 



(28) 



Aix^Zaff^Za) 

where the function fi(Za) is related to the total cross section a^f of the bound-free photo 
production, 



a bf = 4TTa (Za) 3 fx(Za) 



m 



(29) 



The function fi(Za) is of the order of unity for all values of Z, see Ref. jll|]. Since the 
Coulomb correction da@ /de± has opposite sign for electron and positron in the region 
where both particles are relativistic, the correction cr^ is determined by the region where 
momentum of electron or positron is of the order of m. The correction a^ can be obtained 



from our results using the relation 



a 



(i) 
c 



de. 



da c da c , - 8a 

+ i—{Z ->■ -Z) + f{Za) 



de. 



de. 



(30) 
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Our numerical data are in a qualitative agreement with Eq. (|28|) though accuracy is not 
very high because of strong cancellations between all terms in the integrand in Eq. (1301) . 

We conclude this section by discussing how it is possible to apply our results for photo- 
production on heavy atoms, i.e., how the effects of screening can be accounted for. As it was 
pointed out above, the main contribution to the high-energy photoproduction cross section 
close to the end of spectrum of electron (or positron) comes from distances r ~ 1/m r scr , 
where r scr w 1/ [mZ l ^a) is the typical screening radius. Therefore, one can account for the 



screening by employing the prescription formulated in Refs. [23l425l|. Namely, the spectrum 
of slow electron in the screened Coulomb field can be obtain from that in unscreened field 
by means of the shift e_ — > e_ + A and values of the energy shift A for various atoms 



are presented in Refs. 



23 



3M25I]. Analogously, for slow positron the corresponding shift is 
e + — > e + — A. In all cases A/m < 4 x 10 -2 . Thus, the effect of screening in our problem is 
important only for very small electron or positron velocities. 



VI. CONCLUSION 

We have calculated exactly in the parameter Za the cross section of e + e~ photoproduction 
in a Coulomb field at u ^> m and e_ > m (slow electron) or e + > m (slow positron). In the 
wide region, our results differ essentially from those obtained in the Born approximation as 
well as from the results which take into account the Coulomb corrections obtained at e_ ^> m 
and e + ^> m. Therefore, the Coulomb correction to the spectrum can be approximated by 
its high-energy asymptotics only at rather large to {to > 30 m). We have found that the 
cross section of bremsstrahlung in a strong Coulomb field by a high-energy electron in the 
region where a final electron has the energy t\ > m coincides with the cross section of e + e~ 
photoproduction at e_ > m (slow electron). Finally, we have seen that the effect of screening 
for photoproduction close to the end of electron (positron) spectrum is important only for 
very small velocity of electron (positron). 
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Appendix A: Wave functions and Green's function 

The positive-energy wave function L^ )M (p, r) with total angular momentum j, parity 
(— iy +a / 2 (a = ±1), and projection fj, of the total angular momentum on some quantization 



axis has the form 



13| 



Uj,w[p, r) = — 

' ' -o g(r)n j>v>ll (n) 



f(r) = ^jl^e^ ^^ (2pr)^ Im {e^F( 7 - iv, 2 7 + 1, -2ipr)} , 

g(r) = y^ e (-/ 2 ) ^^Tp ( 2 ^) 7 Re {z l(pr+ °ni ~ i", 27 + 1, -2ipr)} , 

, . cr . a Zae /. 1\ /— = 

i=J + -' /= J~-' ^ = — , K = cr|j + -1, 7 = a/k 2 - (Za) 2 , 



2 ' " 2 ' p 

z/fe — m 



7T 

4 = (1 — cr)— + arctan 



a = ±1 , e 4 = — , n — - (Al) 

'y + iv r 



e(7 + k) 

where F(a, fl-,x) is the confluent hypergeometrical function and f^ jjU (n) is a spherical 
spinor. The negative-energy wave function Vj )(7)jU (p, r) employed here can be obtained from 
Uj^nip, r) by the replacement e — > — e. 
If e ^> m and r ~ 1 /m then pr 3> 1 and 

f(r) = sin(pr — In/ 2 + Za ln(2pr) + S K ) , 
g(r) = cos(pr — In/ 2 + ln(2pr) + 8 K ) , 

7T 

<5 K = £ + (Z - 7)- - argT(7 + 1 + iZa) , 

x 7T / Za \ . . . 

£ = 1 - a - + arctan — — . (A2) 
4 V7+ |«|/ 

The high-energy asymptotics of the functions fir) and g(r) for negative-energy states (e <C 
— m) are given by Eqs.( 1A2l) with the replacement Z — > —Z. 

Let us consider the case e — > m, v = p/e <C 1. At |y| — >■ oo and fixed x, 



r(a? + i2/)| ->• v^e" Ws/|/2) \y \ x ~ 1/2 . (A3) 
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Then, at y — > oo, x — > 0, and fixed u = xy we have 

r(2 7 + 1 



F(7 — iy, 27 + 1, — ix) — )■ 



J 2j (2y/u) + — [mJ2 7+2 (2-v/m) - 27 V / mJ2 7+ i(2 1 /u)] 
^2/ 



(A4) 

where J v (x) are ordinary Bessel functions. As a result, we find the following expressions of 
the low-energy asymptotic of the functions /(r) and g(r) 



f(r) = - 7)^r(2Vw) + \/^J 27+ i(2v^)] 



g(r) = aVwrZaJ 27 (2\/u) , w = IZamr . 



(A5) 



It is seen that both /(r) and g(r) are proportional to y/v and thus are of the same or- 
der at Za ~ 1. However, if Za <C 1 then /(r) ^> g(r) at w ~ 1. Using the relation 
Ju+i(x) = iy/x)J u {x) — J' u (x), where J' u (x) = dJ v {x)jdx^ we can write the expression for 
the asymptotics of the function f(r) in the convenient form 



f(r) = <y\j^ [«J 27 (2V^) - ^J'^(2Vu)] 



(A6) 



If e — > —m , v = p/\e\ <C 1, then the functions f(r) and g(r) are exponentially small. We 
have for Za/r ^> mv 2 : 



f(r) = -v^e-^/j^^) 

g(r) = 



V7T 



-(ivZa/v) 



Za 



[(« + 7 )/ 27 (2v^) + v/U/ 27+1 (2Vw)] 



(A7) 



In the Coulomb field, the discontinuity of the quasiclassical electron Green's function on 
the cut reads 

+00 



— T2- exp 
smh r 



i2Za-r + zp(ri + r 2 ) cothr 
P 



x 



7°e + m + -7 ■ (ni — n 2 ) coth r 



Jo(w) 



iJi(w) 



w 



p (r 2 - ri) 
2 sinh 2 r 



Zarwy 



x 7 ■ (ni + n 2 ) - Zap coth r 7 [1 - (7 • n 2 )(7 ■ rij)] 

n,2 



»1,2 



10 



p- v /2rir 2 (l + ni • n 2 ) 



(A8) 



ri j2 sinh r 

In Eq. (1A8j) . the contour of integration over r passes in the positive direction around the 
point t = 0. If p ^> m, ri ~ r 2 ~ 1/m and (1 + n x • n 2 ) ^> m 2 /p 2 , then the argument of the 
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Bessel functions is large and Eq. (1A8j) becomes essentially simpler: 

V 



bG{x,y\\e\) 
8G{x,y\-\e\) 



2ttR 
P 



[ii 3 sm{pR + 2Zas) - 7 • n cos{pR + 2Zas)} 



n 



2nR 
x-y 
R 



[^7° sm(pR — 2Zas) + 7 ■ n cos(pR — 2Zas)] 
x + y + R 



In 



y/2(xy + x-y) 



R 



\x -y\. 



(A9) 
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u- plane 



-1 



1 



FIG. 1: Contour of integration in the complex plane of the variable u used to perform the integral 
inEq. (USD. 
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FIG. 2: The cross section ua~ 1 da/de- of e + e~ pair production at zero electron velocity in units 
a = a(Za) 3 /m 2 . Solid curve: our results via Eq. (f!4j) . dashed curve: the results of Ref. [20] (see 
also Eq. (p~6|) ) . dotted curve: the results of Ref. Q] obtained at uj = 40MeV and /3_ = 0.1265. 
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FIG. 3: The dependence of the Coulomb corrections uio^doc / 'de- (continuous line) with o$ = 
a(Za) 2 /m 2 on the scaled electron energy e_/m at different values of Z. The dashed line represents 
the leading-order Coulomb corrections in the limit e_ 3> m, while the dotted line also includes 
corrections proportional to m/e_. The dashed and the dotted lines start at e_ = 5 m because the 
corresponding asymptotics are valid at e_ S> m. 
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FIG. 4: The dependence of the Coulomb corrections uja^doc / de+ (continuous line) with o$ = 
a(Za) 2 /m 2 on the scaled positron energy e+/m at different values of Z. The dashed line represents 
the leading-order Coulomb corrections in the limit e+ ^> to, while the dotted line also includes 
corrections proportional to m/e + . The dashed and the dotted lines start at e+ = 5 to because the 
corresponding asymptotics are valid at e + S> to. 
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